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ABSTRACT 

Context. High resolution observations of the asteroids Iris and Juno have been performed by means of the UVES spectrograph at the 
ESO VLT to obtain the elfective accuracy of the spectrograph's radial velocity. The knowledge of this quantity has important bearings 
on studies searching for a variability of the fine structure constant carried on with this instrument. 

Aims. Asteroids provide a precise radial velocity reference at the level of 1 m s"' which allows instrumental calibration and the recog- 
nition of small instrumental drifts and calibration systematics. In particular, radial velocity drifts due to non uniform slit illumination 
and slit optical misalignment in the two UVES spectrograph arms can be investigated. 

Methods. The position of the solar spectrum reflected by the asteroids are compared with the solar wavelength positions or with that 
of asteroid observations at other epochs or with the twilight to asses UVES instrumental accuracy . 

Results. Radial velocities offsets in the range « 10-50 m s"' are generally observed likely due to a non uniform slit illumination. 
However, no radial velocity patterns with wavelength are detected and the two UVES arms provide consistent radial velocities. These 
results suggest that the detected Aor/a variability by Levshakov et al. (2007) deduced from a drift of -180 ± 85 m s"' at Zabs =1.84, 
between two sets of Fe ii lines falling in the two UVES arms may be real or induced by other kinds of systematics than those in- 
vestigated here. The proposed technique allows real time quality check of the spectrograph and should be followed for very accurate 
measurements. 
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1. Introduction 

Radial velocity precision is required in several fields of astro- 
nomical research ranging from the detection of exoplanets to the 
study of the variability of the fundamental physical constants. 
To reveal the presence of an orbiting planet dedicated spectro- 
graphs have been manufactured to achieve the best accuracy in 
the radial velocity. With HARPS at the 3.6 m telescope a rela- 
tive precision of 1 m s ' or higher has been achieved when the 
full optical stellar spectrum of a solar type star is recorded and 
compared in different epochs. Search for a possible variability 
of the fine structure constant, Aa/a- (a^ - a)/a, at a redshift 
z, is currently carried out by measuring line shifts between dif- 
ferent lines of absorbers observed in spectra of distant quasars 
which show different sensitivities to a (Webb et al. 1999, Dzuba 
et al. 2002). QSOs are rather faint and require large telescopes 
such as VLT or Keck combined with the high resolution spec- 
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trographs, UVES and HIRES respectively. Murphy et al. (2004) 
claim that Aaja- -5.7 ± 1.1 ppm (ppm stands for parts per mil- 
lion, 10"^) by averaging over 143 absorption systems detected 
in HIRES/Keck telescope spectra of QSOs in a redshift range 
0.2 < z < 4.2 implying that in the past the fine structure constant 
was smaller. On the other hand no variability has been measured 
by a different group at the VLT with UVES adopting similar 
techniques (Quast et al. 2004, Chand et al. 2004, Levshakov et 
al. 2005, 2006; but see also Murphy et al. 2006 and Srianand 
et al. 2007). More recently Levshakov et al. (2007) measured 
a radial velocity difference of -180 ± 85 m s"' between Feu 
transitions falling in the two different arms of the UVES provid- 
ing evidence for a variation in the fine structure constant Aa/a = 
5.4 ± 2.5 ppm, with the fine structure constant being larger in the 
past at odds with what found by Murphy et al. (2004). Given the 
importance of these results for fundamental physics a thorough 
investigation of systematic errors to rule out possible instrumen- 
tal shifts which may occur during UVES observations is rather 
crucial. 

Spectroscopic observations are generally calibrated in wave- 
length by means of standard calibration lamps, namely the ThAr 



2 



P. Molaro et al.: UVES radial velocity accuracy from asteroid observations 



lamps. However, to achieve a Ao'/o' of 1 ppm a precision of 30 
m s"' in the radial velocity of the most sensitive lines is re- 
quired challenging the spectrograph precision. Small instrumen- 
tal effects could be present since the light paths of calibration 
and stellar beams are different when entering the spectrograph 
slits. Instrumental flexures, temperatures and atmospheric pres- 
sure instability can produce small radial velocity shifts between 
calibration and science observations. Temperature and pressure 
variations as small as Ar = 0.3 K or a AP - 1 mbar produce 
a drift of ^50 m s"^ (Kaufer et al. 2004). These effects can be 
minimized by taking ThAr lamps immediately before or after 
the science exposures if ambient conditions do not change in 
the meantime. However, an uneven illumination of the slit may 
cause spectral shifts and therefore errors in the measurements of 
radial velocities. This problem is particularly acute in the case 
of UVES observations with the dichroic mode where the light 
enters two distinct slits of the two arm spectrograph. Possible 
effects of different illumination of the two slits of the blue and 
red arms of UVES are unknown. 

To probe small possible instrumental effects in UVES, we 
observed the solar spectrum reflected by asteroids which are 
sources with radial velocities known at the m s ' level. This 
accuracy is not required by the majority of the observations but 
is crucial for the investigation of variability of the fine struc- 
ture constant. The presence of a variability of fundamental di- 
mensionless constants would be a discovery of the outmost im- 
portance in theoretical physics with far reaching implications 
(Copeland et al. 2006, Avelino et al. 2006, Martins 2006, Fujii 
2007). 



2. Observations and data analysis 

Observations of two asteroids Iris and Juno together with obser- 
vations of the sunlight at twilight have been collected at the VLT 
by means of the UVES spectrograph between December 2006 
and January 2007 as reported in Table [T] UVES is a two-arm 
crossdispersed echelle spectrograph with the possibility to use 
dichroic beam splitters and to record most of the optical spec- 
trum with one observation (Kaufer et al. 2004). The observations 
have been taken with the dichroic mode, the ESO DICl, allow- 
ing simultaneous observations of the blue and red arms. This 
dichroic has a cross-over wavelength at 450 nm and the central 
wavelengths were set at 390 nm for the blue and at 580 nm for 
the red arms respectively, allowing a full spectral range from 350 
nm up to 680 nm. 

We used a 0.5 arcsec slit providing a resolving power of 
about X 80, 000 which is the maximum resolution that can 
be reached with still adequate sampling of the PSF. It is relevant 
to note that the target is centered on one of the spectrograph slits, 
the red slit in our case, while there is no way to check the optical 
centering on the blue arm slit, directly. The slits were aligned 
with the parallactic angle in order not to miss light due to at- 
mospheric diffraction. The two UVES arms are equipped with 
CCD detectors, one single chip in the blue arm and a mosaic 
of two chips in the red arm. The blue CCD is a 2Kx4K, 15 jj. 
pixel size thinned EEV CCD-44 while the red CCD mosaic is 
made of an EEV chip of the same type and a MIT/LL CCID-20 
chip for the redder part of the spectral range. Each arm has two 
crossdisperser gratings working in the first spectral order 

Asteroids are apparently fast moving objects and the geocen- 
tric radial velocity changes by about 1ms"' in about 4 minutes 
typically, thus limiting the maximum exposure time also with the 
largest telescopes. At the epoch Iris and Juno were of 8 and 10.6 



mag, respectively, and we kept the exposures at 300 and 900 s 
achieving a signal-to-noise between 100-200. 

The observations were bracketed by ThAr standard calibra- 
tion lamps. Calibration and science observations were carried 
on in the attached mode to avoid the automatic resetting of the 
spectrograph position, implemented by ESO on 26 Dec 2001, 
to compensate thermal drifts in the dispersion direction between 
daytime calibration frames and science observations. The auto- 
matic resetting of the instrument allows calibration frames to be 
taken in daytime with an economy in terms of observing time, 
but it makes an accurate calibration problematic. 

The data reduction have been performed by means of the 
UVES Pipeline in the MIDAS echelle context. The wavelength 
calibration has been performed using the new atlas of ThAr spec- 
trum by Lovis & Pepe (2007), which increases the laboratory 
wavelength precision by means of HARPS observations, with 
the line selection suggested by Murphy et al. (2007) to avoid 
blends. Mean residuals of < 0.37 mA for the blue arm, < 0.46 
mA for the red low and < 0.55 mA for the red up are gener- 
ally obtained providing a velocity accuracy at the central wave- 
lengths of < 25 m s ' in the red and of < 30 m s ' in the blue as 
shown in detail in Fig. 1. We note that these residuals are about 
one order of magnitude smaller than those derived by Chand et 
al. (2006). They can be further improved for limited portions of 
the spectrum where the reduction is optimized as achieved by 
Levshakov et al. (2007). As shown by de Cuyper & Hensberge 
(1998) an accuracy of 10"^ of the pixel, which in UVES is a!l5 
m s ' , is attainable for non blended ThAr lines with more than 
10^ detected electrons in the central pixels. However, the accu- 
racy of the ThAr lines themselves is of the order of 10 to 100 
m s ' and this is not directly reflected in the residuals of the 
wavelength calibration. The reduced spectra have been normal- 
ized manually tracing the continuum by means of the standard 
MIDAS routine. 

For the reduction of the twilight spectra we skipped the au- 
tomatic sky subtraction and extracted the spectra from the cal- 
ibrated frames manually by using the standard MIDAS echelle 
commands, using a slit height of 8 pixels to minimize effects 
due to the small curvature of the slit projection on the detector. 
To check the curvature effects we extracted 2 spectra from the 
same image with an extraction slit of 2 pixels and offset by 3 
pixels above and below the central position of the order. The po- 
sition of spectral lines on the two extracted spectra did not reveal 
notable shifts due to curvature effects. 

The angular sizes of the two asteroids in the epoch of ob- 
servations were of 0.278 arcsec and 0.263 arcsec, and always 
smaller than the night seeing. They are effective point sources 
and the light follows the same path through the atmosphere, 
telescope and spectrograph not differently from a QSO or other 
point-like sources. Thus with asteroid observations the radial ve- 
locity accuracy could be monitored along the echelle orders for 
the whole frame in a much better way than with the calibration 
lamp since the asteroid lightpath takes into account the atmo- 
spheric variations and the centering of the object on the slit. In 
particular, in the case of UVES observations which make use 
of the dichroic we can monitor the response of the two separate 
arms. In the twilight spectrum the diffuse day-light illuminates 
the slit uniformly so that a comparison between the radial ve- 
locity of the asteroid and the day-light probes slit illumination 
effects on radial velocities. 

UVES is also equipped with an iodine absorption cell which 
can be inserted in the beam to obtain a dense grid of iodine 
absorption lines superimposed on the target spectrum. The io- 
dine cell currently mounted on UVES produces a rich absorp- 
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Table 1. Journal of asteroid observations and basic data. Ceres spectra has been observed by HARPS. Expected radial velocities and its compo- 
nents are given in columns 6, 7 and 8. For Ceres the values refer to mid exposure. 



Name 


Date 


JD 


V 


exp. 


R^ast-par 




ARV 








mag 


sec 


km s ' 


km s ' 


km s"' 


Iris 


18/12/06 


2454088.518086 


7.83 


300 


12.707 


1.704 


14.411 




22/12/06 


2454092.514539 


7.95 


300 


13.777 


1.830 


15.607 




23/12/06 


2454093.514000 


7.98 


300 


14.030 


1.862 


15.892 




24/12/06 


2454094.515400 


8.01 


300 


14.283 


1.893 


16.176 




25/12/06 


2454095.513566 


8.04 


450 


14.521 


1.924 


16.445 


Juno 


24/01/07 


2454125.880451 


10.62 


600 


-21.403 


3.731 


-17.672 




25/01/07 


2454126.873002 


10.61 


900 


-21.364 


3.721 


-17.643 




29/01/07 


2454130.873891 


10.57 


900 


-21.065 


3.683 


-17.382 




31/01/07 


2454132.846757 


10.55 


900 


-20.949 


3.664 


-17.285 


Ceres 


15/07/06 


2453932.837256 


8.00 


1800 


-11.288 


0.456 


-10.832 




22/05/06 


2453877.919808 


8.85 


900 


-22.707 


0.690 


-22.017 



Table 2. Sky observations and data, * spectra taken with HARPS 



Date 


JD 








km s"' 


18/12/06 


2454088.480700 


0.238 


22/12/06 


2454092.481400 


0.263 


23/12/06 


2454093.479390 


0.272 


24/12/06 


2454094.478500 


0.279 


25/12/06 


2454095.480300 


0.286 


25/01/07 


2454126.472293 


0.573 


31/01/07 


2454132.482003 


0.629 


14/07/06* 


2453931.425152 


0.302 


22/10/05* 


2453666.426383 


-0.070 



3. Asteroids as radial velocity standards 




Fig. 1. Typical order residuals of the wavelength calibration, 
namely the difference between the measured and laboratory 
wavelength of the ThAr lines used in the calibration. The plotted 
ones are for Iris 22 Dec 2006. The three groups, with mean val- 
ues and 1 cr dispersion over plotted, refer to the 3 CCDs which 
have been reduced independently. 



tion line spectrum in the range 490-640 nm. Butler et al. (2004) 
achieved an accuracy of 0.42 m s"' for UVES with observations 
of a Cen A by means of a series of 3013 spectra of 1-3 s expo- 
sures, but after correcting for trends and jumps. However, iodine 
cell is not well suited for measuring accurate positions of QSO 
absorption lines which fall very far apart, and we are not aware 
of its use for this purpose. 



Out of 1 1 1 stars observed in 20 years with the two CORAVEL 
spectrometers only a minor fraction shows variability of a!200 
m s"' (Udry et al. 1999). Thus radial velocity standard stars 
provide a reference system of radial velocities with a precision 
of several hundred m s"' . Among the celestial sources the aster- 
oids are probably the best radial velocity standard sources and 
at least for two reasons. The former is that they reflect sunlight 
without any modification of the solar spectrum and the latter is 
that their velocity component with respect to the observer can 
be predicted with very high accuracy reaching the m s"' level 
(Zwitter et al. 2007). 

The first condition is strictly valid only for relatively large as- 
teroids with a nearly spherical shape which produce a constant 
reflectance of the sunlight. The two selected asteroids Iris and 
Juno have radii of 99.9 and 1 17.0 km respectively and a spheri- 
cal shape. On the 18 Dec observation of Iris the illuminated frac- 
tion were of 97.07% and on the 24 Jan 2007 Juno had a 97.26 % 
reflectance so that the reflected and the direct solar spectra are 
likely identical. Variation of reflectance with wavelength or pres- 
ence of regolith developed by meteoroid impact on the asteroid 
do not affect high resolution spectra. Also the asteroid rotation 
does not affect the solar spectrum and is much smaller than the 
solar one. The rotational periods for Iris and Juno are of 7. 14 and 
7.21 hours respectively. Thus their rotational velocity would be 
of about 25 m s"' , which is much lower then the solar one and 
it will not cause further significant broadening. 

The second reason is that the component of their motion rel- 
ative to the observer on the earth can be calculated with extreme 
accuracy. For asteroids with radar monitoring, the orbital com- 
putations can take into account the interferences of other bodies 
of the solar system including the major asteroids and reach pre- 
cisions at the level of the m s"' (Zwitter et al. 2007). 
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Table 1 reports the motion components and the resulting ex- 
pected radial velocity shifts ARV. Ephemeris for our objects has 
been computed by means of the JPL's Horizons systenr] which 
provides accurate ephemeris for the minor bodies of the solar 
system. The sunlight reflected by the asteroid is shifted by the 
heliocentric radial velocity of the asteroid with respect to the 
sun at the time ti when the photons left the asteroid and were 
shifted by the component of the earth rotation towards the aster- 
oid at the time t2, when the photons reach the earth. The latter 
is the projection along the line-of-sight of the asteroid motion 
with respect to the observer at the Paranal site adjusted for aber- 
ration, and comprises both the radial velocity of the asteroid and 
the component due to the earth rotation towards the line of sight. 
At Paranal the observed asteroid radial velocity is 



16.5 



ARV ^ (RVas,-par + RV„„^q) 



(1) 



We also take as reference several twilight spectra which are 
listed in Table 2. The radial velocity of the skylight reflected by 
the terrestrial atmosphere is also shifted by the heliocentric earth 
radial velocity. This component, RVp„i._Q, is given in the last col- 
umn of Table 2. The precise position of the scattering of solar 
light by the atmosphere is not known but it should be within 10 
km from the ground. In the next sections we will show that the 
scattered light from the atmosphere probably keeps the transver- 
sal motions of the atmosphere and therefore it is not possible to 
predict its velocity with the desired accuracy. 

4. Asteroids with UVES 

4.1. Solar absolute reference 

The highest quality solar spectra in the optical domain are the 
FTS solar flux and disc-center atlas obtained at the McMath 
telescope at Kitt Peak by Kurucz et al. (1984) and Brault & 
Neckel (1987). These atlases achieve a signal-to-noise ratio of 
about 2500 with a resolving power of 400,000. AUende Prieto & 
Garcia Lopez (1998a,b) used these atlases to measure the cen- 
tral wavelength for a considerable number of lines. Gravitational 
shifts and convective motions are responsible of line to line dis- 
placements which can be of several hundreds of m s ' . These 
displacements vary with the solar cycle showing a modulation 
with a peak to peak variation of 30 m s"' on the 1 1 years solar 
activity period with the positions more redshifted in correspon- 
dence of the maximum of activity (Deming & Plymate 1994). 
However, McMillan et al. (1993) did not reveal any drift within 4 
m s ' in the solar line position from a long data series spanning 
the period from 1987 to 1992. Allende Prieto & Garcia Lopez 
(1998a,b) line positions have a precision of the order of a!50-150 
m s ' so that they provide absolute reference at this level. The 
lines formed at the top of the photosphere show shifts close to the 
gravitational redshift of 636 m s ' while the other lines show the 
effects of convective motions with variable blue shifts of several 
hundreds of m s ' . Lines with equivalent width stronger than 
200 mA are rather insensitive to the convective shifts and have 
been used to estimate the absolute zero of the scale. The value 
at the plateau level is of 612 ± 58 m s"' in the case of the solar 
atlas of Kurucz et al. (1984) which shows the results closest to 
the theoretical gravitational shift. 

We thus compare the measured line positions of the aster- 
oid spectra with the solar line positions provided by Allende 
Prieto & Garcia Lopez (1998a) for the solar atlas of Kurucz et 
al. (1984). In fact the sun light reflected by the asteroids is a 
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4000 4100 4200 4300 4400 4500 

X(A) 

Fig. 2. Line shifts of Iris 23 Dec 2006 with reference of the 
Allende Prieto & Garcia Lopez (1998a) solar line wavelengths, 
see text for details. The dotted line shows the expected veloc- 
ity of the asteroid. The top panel refer to the red-up CCD, the 
middle panel to the red-low CCD, and the bottom panel the blue 
CCD. The mean values and their dispersion are shown with the 
squares in the middle of each panel. 



sort of integrated solar flux as the Kurucz et al. (1984) atlases. 
Fig. 2 shows the ARV measures for the Iris spectrum of 23 Dec 
2006. The figure shows that there are no major wavelength cali- 
bration inaccuracies at the level of 200 m s"' which corresponds 
to about 0.1 of the pixels size. The result shows a mean value of 
ARV - 15.614 + 0.203 km s"' for the 75 fines measured in tfie 
red-up CCD, a ARV = 15.668 + 0.234 km s ' for 96 lines in tfie 
red-low CCD, and of 15.582 + 0.300 km s ' for 63 fines in tfie 
blue CCD. Considering tfiat tfie expected velocity is of ARV = 
15.620 km s ' , tfiere is an excellent agreement witfi tfie red-up 
CCD and a sligfit offset of about 50 m s"' and 30 m s"' witfi 
tfie red-low and blue CCD respectively. Despite tfie scatter of 
~200-300 m s"' , tfiis analysis sfiows tfiat tfiere is no significant 
offset between tfie two arms of UVES implying tfiat tfiere is no 
mis-centroiding of tfie target on tfie two slits of UVES arms. 

4.2. Asteroid versus asteroid 

To overcome tfie intrinsic uncertainties in tfie positions of tfie 
solar lines we compare tfie solar spectrum from an asteroid taken 
in two different epocfis. In tfiis way eacfi line is compared witfi 
itself leaving only tfie instrumental and calibration imprinting 
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on the change of the asteroid radial velocity between the two 
epochs. 

To measure accurately the radial velocity difference between 
two lines we have implemented and adapted a procedure from 
Levshakov et al. (2006). The most probable ARV between two 
lines is found by varying ARV by small incremental steps and 
estimating the of the fit. Fig. 3 shows a portion of the red- 
low frame of Juno 31 Jan and of the sky spectrum of the same 
day around line 5217.09 A. The S/N are of 126 and 295 for Juno 
and twihght, respectively calculated from two nearby continuum 
windows bracketing the line position. The range used in the fit- 
ting is marked by thick curves on the upper panel of the figure 
and we consider only the central parts of the absorption fines to 
avoid the influence of the wings. 

From the fit of points in the vicinity of the global x^^^ '■^^ 
procedure computes a parabola with the radial velocity differ- 
ence as variable. The Icr uncertainty interval is then calculated 
from the parabola when;^f^(Av) - Xmin ~ ^- particular 
case we have obtained ARV = -17.723 ± 0.033 km s~^ at Icr. 

The error is rather typical of our measurements and corre- 
sponds to about 0.02 of the pixel size. For instance, in comput- 
ing the difference between the Iris spectra taken on 18 and 22 
Dec 2006, the 203 lines measured have a mean error of 39 ± 3 
m s"' , and of 37 + 2 m s"' , respectively. The error is mainly 
the photon noise error and it depends from the signal to noise 
ratio of the two spectra. This error sets the precision of our anal- 
ysis and the level of instrumental effects which can be recovered. 
In principle with higher signal-to-noise spectra this level can be 
further improved. Wavelength calibration errors are not expected 
to contribute very much to this error because even if the wave- 
length cafibrations of two spectra are performed independently, 
they likely make use of the same Ar or Th lines in deriving the 
calibration coefficients. 

The results of the radial velocity difference between the Iris 
spectra taken on 18 with those of 22 and 23 Dec 2006 are shown 
in Figs. 4 and 5. The measures are performed on lines falling 
onto 7 orders for each CCD frame selected to map the full spec- 
tral range. On the right side of each panel the average value 
for each single order is reported with the sample standard de- 
viation. As it can be seen from the top panel of the figure, the 
measurements do not show evidence for trends within an indi- 
vidual order, and the measures are normally distributed around 
their mean value. At the bottom of the figure the mean values for 
each order are plotted as a function of the order number. There 
is no evidence of any pattem of the measured radial velocity 
with wavelength from 3500 A up 6750 A with measurements 
involving 3 CCDs and two spectrograph arms. For the 18-22 
comparison the mean of the 3 CCDs are ARV= 1.157 ± 0.073 
km s"' for the blue, 1.140 ± 0.048 km s"' for the red-low and 
1.134 + 0.048 km s"' for the red-up. The excess in the dispersion 
observed within each order reflects the combined contribution of 
wavelength calibration and data reduction errors with the statis- 
tical error. The mean of the mosaic of the two CCDs of the red 
arm is only 19 ms"' away from the value of the blue arm. For 
the comparison between 18 and 23 Dec observations we have a 
mean value for the blue CCD of 1.487 ± 0.042 km s'^ and for 
the mean of the two red CCDs a value of 1.464 ± 0.036 km s"^ , 
or 23 m s"' away from the blue arm. 

Therefore, there is no evidence for a signiflcant mis- 
aligrmient between the two arms of the UVES spectrograph. 
However, the expected radial velocity difference between the 
two epochs in which the asteroid was observed is of 1.190 
km s ' and of 1. 447 km s ' , respectively. Taking the mean 



5 o.a 



"date: Jan 31, 2007 
















: : : , : : 1 1 






^ --- Juno 

Sky 



Wavelength (A) 




17.69 17.7 17.71 17. 7^ 17 73 17,74 17.75 
Av, km s"' 



Fig. 3. Procedure for the determination of an accurate ARV and 
its error. On the top panel the normalized sky spectrum is shghtly 
shifted vertically for display purposes. 

value of the two arms we miss the expected velocity by 43 and 
24 m s"' in the two comparisons. This implies a sort of sys- 
tematic error in one or in all the observations. The origin of this 
systematic error is likely to be ascribed to a non-uniform iUumi- 
nation of the slit. 

Given that there is no evidence for a systematic behavior 
within the orders in the rest of our observations we have per- 
formed a cross-correlation to get order shifts by means of the 
IRAF-rvsao XCSAO routine. For this kind of analysis particu- 
lar care has been adopted in selecting spectral regions without 
telluric lines which perturb the cross correlation. At the bottom 
of Figs. 4 and 5 the measures based on single fines, plotted in 
dots, with those performed by means of the XCSAO, plotted 
in diamonds, are showing that the two procedures are provid- 
ing consistent results. The results are reported in Table 3 and 
summarized in Table 4. These measurements performed on the 
whole set of observations at our disposal confirm that there are 
not notable patterns with wavelength, no offsets between the two 
UVES arms and that offsets with respect to the expected velocity 
in the range 10-50 m s ' are common. 

5. Asteroids with HARPS 

To check the whole procedure by means of a different instru- 
ment specifically designed for high precision radial velocity 
studies, we retrieved two reduced spectra of Ceres from the pub- 
lic HARPS archive and applied the same kind of measures per- 
formed with UVES. The High Accuracy Radial velocity Planet 
Searcher at the ESO La SUla 3.6m telescope is a spectrograph 
dedicated to the discovery of extrasolar planets through radial 
velocity oscillations. It is a fibre-fed high resolution echelle 
spectrograph and is contained in a vacuum vessel to avoid spec- 
tral drift due to temperature and air pressure variations. There 
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Fig. 4. Radial velocity difference from Iris 18 and 22 Dec 2006. Residuals correspond to the difference Iris(22) - Iris(18). The 
predicted ARV is 1.190 km s"^ . Individual echelle orders (numbered by bold) are shown in the upper panels. For each order the 

mean value Av and the sample standard deviation cr are indicated. These values Av and cr are also shown by dots with error bars in 
the corresponding low panels where results obtained through the cross-correlation analysis (diamonds) are plotted for comparison. 



are two fibers, one collects the star light, while the second is 
used to record simultaneously a ThAr reference spectrum. Both 
fibres are equipped with an image scrambler to provide a uni- 
form spectrograph pupil illumination, independent of pointing 
decentering. In this way the instrument is able to obtain a long 
term radial velocity accuracy of the order of 1 m s~^ for the 
entire optical spectrum of a slow rotating G type star or cooler 
(Pepe et al. 2005). HARPS has a resolving power of R « 120,000 
and provides a sampling of the slit of FWHM = 4.1 pixels of 
15 fi size. Due to the relatively smaller size of the telescope the 
exposures are rather long, being namely of 1800 s and 900 s 
(see Table 1). In the course of the exposure the radial velocity 
of the asteroid changes by =^50 and 25 m s ' respectively. The 
expected velocities reported in Table 1 refer to the mid exposure 
times. 



In Fig. 6 the radial velocities measured between the observa- 
tions of Ceres taken on 22 May 2006 and on 15 July 2006 are 
given. The accuracy of the measure of a pair shift is now better 
than a!20 m s"' and the line to line variation of the positions is 
almost entirely due to errors in the wavelength calibration. The 
mean of the blue CCD Linda is ARV ^ 3.072 ± 0.010 m s ' , 
and the mean of the red CCD Jasmin is 3.080 ± 0.010 m s"' . 
The predicted ARV shift is of 3.070 km s~^ and is found in ex- 
cellent agreement with the measured velocity within few m s"' . 
This is highly suggestive that the systematic offset observed in 
the UVES spectra is related to the slit acquisition mode which 
remains the most significant observational and technical differ- 
ence between the two spectrographs. 
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Fig. 5. Same as Fig. 4 but for the comparison between Iris 18 and 23 Dec 2006. The predicted ARV is 1.447 km s ^ Residuals 
correspond to the difference Iris(23) - Iris(18). 



6. Twilight 

Sky observations differ from point source observations mainly 
because the slit is uniformly illuminated. Thus in principle a dif- 
ferential measure of a point-hke source as an asteroid with the 
same feature from the sky spectrum allows us to probe radial 
velocity drifts induced by a non uniform slit illumination. The 
results of the measures for the Juno observations of 31 Jan are 
shown in Fig. 7. The observations of Juno on 3 1 Jan when com- 
pared with the skylight on 31 Jan show that the A/? V inferred 
from the three CCDs are all consistent with each other. The blue 
CCD gives a mean value of ARV= -17.797 + 0.064 km s ' , 
the red-low a ARV ^ -17.757 ± 0.050 km s~' , and the red-up 
a ARV^ -17.772 ± 0.057 km s"' . The expected velocity is of 
ARV = - 17.914 km s"' , therefore we observed an offset of »150 
m s~^ . This offset is rather high and about a factor three higher 
than that observed in the series of asteroid-asteroid comparison. 
To check the procedure we performed two separate tests. In one 



test we compare the accurate HARPS observations of Ceres with 
a sky spectrum taken with the same instrument, and in a second 
test we compare two twiUght spectra taken with UVES in two 
different epochs. 

Fig. 8 shows the comparison between the spectrum of Ceres 
taken on 15 July 2006 with a sky light taken with the same instru- 
ment on 22 Oct 2005. The mean ARV values in the two HARPS 
CCDs are ARV = 23.657 km s"' for the blue CCD Linda, and 
ARV= 23.647 km s~^ for the red CCD Jasmin, while the pre- 
dicted one is of ARV = 23.709 km s"^ computed for the middle 
Ceres's exposure. This measure shows an offset of a!50 m s~' 
between the Ceres and sky spectrum. This offset is not observed 
when the Ceres of two epochs are compared with each other, as 
we discussed in the previous section. HARPS is fed by fiber op- 
tics and therefore no difference is expected between the two kind 
of measures suggesting that the sky spectrum holds a component 
of motion of several tens of m s ' . This implies that the twilight 
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Fig. 6. A/?y between Ceres 22 May and 15 July 2006 . The ex- 
pected ARVis 3.070 km s'^ 



Fig. 8. ARV between Ceres 15 July 2006 and the Sky on 22 Oct 
2005 . The predicted ARV is 23.709 km s'^ . 



solar spectrum is not a good reference for the determination of 
the zero scale. 

As a second test we compared the sky spectra with each 
other. The difference between the sky spectra taken with HARPS 
on 22 May 2005 and 14 July 2006 are shown in Fig, 9, The ex- 
pected velocity difference for this pair is of -371 m s"' while 
the mean value is -275 ± 39 m s~' . Thus, also in this case we fail 
to reproduce the expected velocity confirming that the sky spec- 
trum is sensitive to unpredictable motions likely due to currents 
in the upper terrestrial atmosphere. 

We also emphasize that close inspection of UVES twilight 
and asteroid solar spectra show that they are not completely 
identical. Small diff'erences at the level of 1-2% are found be- 
tween the twilight spectrum and the asteroid reflected solar spec- 
trum consistently with what found by Zwitter et al. (2007). An 
example of the two spectra, with the skylight Unes shallower, is 
shown in Fig, 10. Similar differences have been found also by 
Gray et al, (2000) and also depending on the angular separation 
from the Sun, According to Gray et al, (2000) the skylight vari- 
ations can be explained as a combination of Rayleigh-Brillouin 
scattering with a second term of aerosol. The measure of the 
FWHM for a representative sample of lines from asteroid spec- 
tra and for twilight are shown in Fig. 1 1 . The FWHM of the 
twilight are significandy larger by about 5 mA with comparison 
to the asteroid fines. The broader twilight fines suggest the pres- 
ence of turbulence in the atmospheric motions which reflect the 
sun light. At twilight a transverse motion in the atmosphere has a 
considerable component in the direction of the sun which is low 





5500 
lambda (A) 



Fig, 9. ARV between HARPS sky of 22 May 2005 and 14 July 

2006, The dotted line is the mean value while the red line shows 
the expected ARV at -371 m s ' . 



above the horizon when the observations were made and produce 
a radial velocity drift when observed in the reflected spectrum. 
A detailed investigation of these effects is beyond the scope of 
this paper, but the presence of this effects shows that the twilight 
spectrum is not a good zero reference point at the level of «100 
m s"^ . 
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Fig. 7. Juno 31 January relatively to the twilight of the same date. Residuals correspond to the difference Juno - Sky. 



7. Implications for Aa/a 

In the Many Multiplet method the measurabiUty of Aa/orfrom 
observations of absorption lines in QSO spectra is based on the 
fact that the energy of each line transition shows a different sen- 
sitivity on a change of a (Webb et al. 1999). Thus, the value 
of Aa/a depends on the measure of the relative radial veloc- 
ity shifts, ARV , between lines with different sensitivity coeffi- 
cients. The relation between the radial velocities and Aa/a is 
(Levshakov et al. 2006): 

Aa 

(V2-Vi) = 2c(Qi-Q2)— , (2) 

a 

where Q is the sensitivity coefficient Q - q/coo, with ojq be- 
ing ffie frequency and q the theoretical so-called q-factor. The 
^-factors have been computed for the most important UV reso- 
nance transitions by Dzuba et al. (2002) and for Fen were re- 
calculated by Porsev et al. (2007). 

The largest A<3 is presently provided by the Fe ii resonance 
lines. By comparing ffie Feii/11608 wiffi and Feii /i2382 or 



^2600 Unes (2x608 = -0.0166, ^2382 = 0.0369, and 22600 = 
0.0367 from Porsev et al. 2007) we obtain lAgl ^ 0.053 which 
is almost two times larger ffian that obtained from a combination 
of offier transitions. In ffiis case a shift of a!30 m s~^ between the 
Fell lines corresponds to a Aa/a of ppm. 

Levshakov et al. (2007) analyzed Feii profiles associated 
with the Zabs - 1 -84 Damped Lya system from UVES observa- 
tions of the quasar Q 1 101-264. The data represent one of very 
few spectra of QSOs obtained with spectral resolution FWHM 
of 3.8 km s ' and S/N > 100. In this work a shift of the relative 
radial velocity between the /1 1608 and /1/12382, 2600 lines of 
ARV= -180 ± 85 m s"' was obtained. With ffie updated sen- 
sitivity coefficients from Porsev et al. (2007) this shift in the 
radial velocity between the Fe n lines corresponds to a Aa/a = 
5.66 + 2.67 ppm. 

The Fen lines fall at i ~ 4566 A and .1 ~ 6765,7384 A, 
respectively, quite far apart in the two different UVES arms so 
that a hidden systematic effect would challenge the interpretation 
as due to variation of a. Levshakov et al. (2007) measured the 
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Table 3. Cross correlation analysis: differential radial velocity shifts measured in km s ' with respect of Iris 18 Dec 2006. The first column reports 
the echelle orders for the three chips blue, red-low and red-up. Xcsao cr are also indicated. 
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31.658 


0.035 


99 


6131.5 


-6193.7 


32.112 


0.045 


32.020 


0.041 


31.795 


0.038 


31.630 


0.038 


100 


6077.0 


-6131.5 


32.197 


0.063 


32.051 


0.059 


31.824 


0.052 


31.589 


0.050 


102 


5995.0 


- 6010.7 


32.184 


0.049 


31.978 


0.042 


31.741 


0.048 


31.610 


0.037 


104 


5841.0 


- 5885.0 


32.172 


0.064 


32.087 


0.064 


31.721 


0.063 


31.601 


0.063 



same velocity between the /12382 and /12600 which is what ex- 
pected since the Q values for A23S2 and /12600 are about equal. 
However, there is no direct way to check out systematic differ- 
ences of the Fen /11608 with lines that fall in the other arm of 
the spectrograph. Different velocity offsets may occur in the blue 



and red frames causing an artificial Doppler shift between the 
Fe 11/11608 and /1/12382, 2600 lines and mimicking a change in 
Aa/a . The set of measures carried out here show that there are 
no ARV offsets between the two UVES arms greater than 30 
m s"^ . This excludes this kind of systematics as a possible ori- 
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Fig. 10. Portion of the asteroid (thick line histogram) and twi- 
light (thin line histogram) spectrum around the line 5217.3 A. 




0.1 0.15 0.2 0.25 

FWHM (A) 



Fig. 11. FWHM for asteroid Juno (solid line) and twilight 
(dashed line). 



gin of the signal detected by Levshakov et al. (2007). Therefore, 
either the detection is real or it is induced by a systematics of 
different kind. 



8. Conclusions 

Observations of asteroids have been conducted with the UVES 
spectrograph at the VLT to probe the radial velocity accuracy 
achievable with the spectrograph. By means of HARPS observa- 
tions we have shown that the asteroid observations are excellent 
radial velocity standards able to probe the instrumental accuracy 
in any particular position of the spectrum down to the limit pro- 
vided by the ThAr wavelength calibration, or 10ms' . 

By comparing the asteroid line positions with the absolute 
ones from solar positions which account for solar convective 
shifts we have shown that the UVES spectrograph is not affected 
by any systematics along the whole optical domain at the level 
where the solar line positions are known, namely of few hun- 
dreds of m s ' . We have further refined the analysis by compar- 
ing of asteroid in different epochs. No major distortions in the 
wavelength are found, namely not higher than about 30 m s"' , 
where this limit is set by the photon noise of our observations. 
We do indeed reveal zero offsets in the range up to ^50 m s ' . 
With reference to similar observations performed with HARPS 
we suggest that this is likely due to a non uniform slit illumina- 
tion. Attempts to use the twilight spectra to quantify the drifts 
induced by non-uniform illumination shows instead that twilight 



spectrum contains additional turbulence and motions, and there- 
fore cannot be used as a reliable zero reference point. 

The recorded spectrum does not show evidence of stretching 
of the wavelength scale or other instrumental effects in excess 
of the uncertainties induced by the wavelength calibration ac- 
curacy. In particular, the two UVES arms which are fed by two 
independent slits do not show signature for radial velocity offsets 
within the present accuracy of 30 m s"' . 

This result has important implications on the search for 
Aa/a currently performed with UVES which relies on relative 
shifts of absorbing lines falling on rather distant spectral regions 
and sometimes belonging to different arms of the spectrograph. 
For instance, Levshakov et al. (2007) measure of a ARV differ- 
ence of -180 + 85 m s"' between Fe ii transitions falling in the 
two different arms of UVES which provides evidence for a vari- 
ation in the fine structure constant Aa/a= 5.66 ± 2.67 ppm. The 
present analysis shows that the line shift is unlikely produced by 
a misalignment of the the two slits at the entrance of the two 
UVES arms. 

The proposed technique has a general validity and can be 
applied to any spectrograph to perform a real time quality con- 
trol of the spectrograph performance during night time while the 
observations are carried on. 
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Table 4. Summary of mean radial velocity shifts measured in km s ' with respect of Iris 18 Dec 2006. 



Iris 


Date 


ARV 


Blue 


cr 


Red-low 


cr 


Red-up 


cr 


Red 


Tot 


cr 


18-22/12 
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0.041 
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0.011 
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-2.055 
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0.015 
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0.042 


Juno 


Date 


ARV 


Blue 


cr 


Red-low 


cr 


Red-up 


cr 


Red 


Tot 


cr 
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32.204 


0.051 


32.101 


0.050 


32.160 


0.041 


32.130 


32.152 


0.063 


18-25/01 


32.054 


32.053 


0.021 


31.983 


0.087 


32.020 


0.052 


32.002 


32.017 


0.063 


18-29/01 


31.793 


31.822 


0.033 


31.730 


0.052 


31.765 


0.047 


31.747 


31.770 


0.058 


18-31/01 


31.696 
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31.590 
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31.579 


31.606 


0.068 
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